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A strong-coupling expansion for models of correlated electrons in any dimension is presented. The 
method is applied to the Hubbard model in d dimensions and compared with numerical results 
in d = 1. Third order expansion of the Green function suffices to exhibit both the Mott metal- 
insulator transition and a low-temperature regime where antiferromagnetic correlations are strong. 
It is predicted that some of the weak photoemission signals observed in one-dimensional systems such 
as SrCu02 should become stronger as temperature increases away from the spin-charge separated 
state. 
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Organic conductors, cuprate ladder compounds and 
High- Tc superconductors are but a few of the condensed 
matter systems currently driving the intense experimen- 
tal and theoretical efforts on strongly correlated electrons 
in low dimension {d — 1 oi d — 2). Angle-resolved 
photoemission experiments (ARPES) on 2D and ID ma- 
terials are beginning to probe the spectral weight 
v4(k, uj), but have not yet given definitive answers to ques- 
tions of prior interest such as spin-charge separation or 
the relation between the Mott transition and antiferro- 
magnetic (AF) correlations. On the theoretical side, the 
Hubbard model (HM) |^ is the simplest one that includes 
the interplay between the strong screened Coulomb re- 
pulsion and the kinetic band energy. In one dimension, 
the exact solution of the HM Q does not allow actual 
calculation of correlation functions, but several numeri- 
cal studies [||-|| have investigated the one-particle spec- 
tral function that is observed in ARPES. Other meth- 
ods (bosonization [p|-pT|, renormalization group Jl2| , p^ , 
conformal field theory |1^) have led to important non- 
perturbative results, but in addition to being restricted 
to the lowest energy excitations, they involve parameters 
that are absent from the microscopic Hamiltonian. For 
the d = cx) case, all quantities of interest can be 

calculated in an essentially exact way, but the extrapo- 
lation to low dimension is problematic. Several attempts 
to develop systematic strong-coupling expansions [ p^j20| 
did not yield the spectral weight. 

The purpose of this letter is two-fold. First, we con- 
struct a strong-coupling perturbation theory that can be 
applied to a number of models in any dimension, and, 
second, we use it to compute the Green function of the 
half-filled HM. This allows us to discuss the Mott tran- 
sition from the viewpoint of the density of states. The 
effects of antiferromagnetic correlations on A(k, uj) are 
discussed, for simplicity, only in ID. Despite the absence 
of phase transitions in ID, the qualitative behaviour of 
yl(k, uj) allows us to define crossovers between regions of 
parameter space where the system behaves somewhat like 
a metal, an insulator or a short-range antiferromagnet. 
The results are summarized by the crossover diagram of 
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FIG. 1. Crossover diagram of the half-filled ID Hubbard 
model with Coulomb repulsion U — 2u and hopping t. 



Fig. 1^, which shows some analogies with one of the pub- 
lished d = OO phase diagrams [0. We conclude with a 
prediction for ID systems of current experimental inter- 
est p. 

First, let us present the strong-coupling expansion it- 
self. Consider a Hamiltonian H = Ji^+Ti} , where the un- 
perturbed part TiP is diagonal in a certain variable i (say 
a site variable), and let us denote collectively by a (say 
a spin variable) all the other variables of the problem. 
This Hamiltonian involves fermions, and is supposed to 
be normal ordered in terms of the annihilation and cre- 
ation operators c^'J . may be written as a sum over i 
of on-site Hamiltonians involving only the operators c^P 
at site i: H° — J^i^ii^la^'^icr) ■ ^or a strong-coupling 
expansion of the HM, H'^ is the atomic limit, namely 
hi[c\^, da) — Ucl-^cl^CiiCii (we will often use u ^ U/2 for 
convenience). We suppose that the perturbation is a 
one-body operator of the form = ^ij^l^Cja ■ 

For the HM, Ti,^ is the kinetic term. Introducing the 
Grassmann field 7icr(''"), 7i*cr (''')) the partition function at 
some temperature T = 1//3 may be written in the Feyn- 
man path-integral formalism: 



Z = J [dj*d-f] exp -jyri^Y. (^) 
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We use the letters (a, b, ...) to denote sets such as (z, a, r), 
for instance : 

•^0 y"(T ah 

A first difficulty arises: There is no Wick theorem 
because Ti^ is quartic instead of quadratic. We solve 
this problem by means of a Grassmannian Hubbard- 
Stratonovich transformation, [2l| ] which consists in ex- 
pressing the perturbation part of the action in Eq. (|^) 
as a Gaussian integral over an auxiliary Grassmann field 
ipia- (r) , ip*^ (t) . Then, the integral over the original vari- 
ables can be performed and Z can be rewritten in the 
form: 

J [drd^p] exp - i^Soir, i'] + E ^i"t[^*' ^] I • 

The action has a free part 

ab 

and an infinite number of interaction terms 



{ai,bi} 



bi..bft 



where the G?f..a« = {lai--lanlt:^-^Oo c 
bi..bR 

nected correlation functions of the unperturbed system. 
The primed summation reminds us that the fields in each 
term share the same value of the site index. We may now 
use Wick's theorem and usual perturbation theory for the 
^'s, the free propagator being F, and the vertices being 
the G^'^'s. The number of auxiliary field propagators de- 
termines the order in y = t for the HM) of a given 
diagram. Finally, the relation between the Green func- 
tion Qab = ~(7a7b) of the original fermions and that of 
the auxiliary field Vab = -~{'>Pa'4'b) ^ matrix form) 
g ^ + V-^VV-^ . If r denotes the self-energy of 
the tp's, one has G = (p-i - V)'^ . 

The above method was applied to the HM 

n = 2uJ24]Al(^^l(^^^ - t (^I<tCjt + H.C.) 

at half- filling up to order t^. The result for Q is a rational 
function of iun'- 



^(k,iw„) 



= 2te(k) 



lUJr, 



u2 {2{iuJnf - U^) 



where d is the dimension of the hypercubic lattice, and 
= Ylm=i cos(fcm). Here we face a second difficulty, 
namely that C/(k, iWn) has pairs of complex conjugate 
poles. This violates the Kramers-Kronig relations and 
leads to negative spectral weight. Note that even in weak- 
coupling theory, truncation of the series for Q leads to 
high-order poles giving negative spectral weight. Since 
we only know Q up to order i^, any function having the 
same Taylor expansion as Q to this order is a priori as 
good an approximation. A physically acceptable solution 
should be causal and have a positive spectral weight, that 
is, be a sum of simple real poles with positive residues. 
We call such a function Lehmann representable (LR). 

In order to obtain a LR approximation, we need the 
following theorem, reported in Ref. |^: A rational iunc- 
tion is LR if and only if it can be written as a finite Jacobi 
continued fraction 
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with bl real and a; > (thereafter conditions CO). 

According to this theorem, the exact Green function 
of any finite system is a Jacobi continued fraction, whose 
coefficients, functions of the hopping t, verify conditions 
CO. If we expand the exact Gj in powers of t to some 
finite order, which is what a strong-coupling expansion 
does, we destroy its continued fraction structure. If in- 
stead we replace ai{t) and bi{t) in Qj by their expansion 
to some finite order, the result should be LR since we ex- 
pect conditions CO to hold for the truncated coefficients 
(at least for t/u small). 

Therefore, to obtain a LR approximation, we seek 
frequency- independent ai{t) and bi{t), such that Qj and 
Q have the same expansion up to order t^. Equating the 
series in t for Q and for Qj at all frequencies determines 
uniquely the leading terms in the t expansion of a/ (t) and 
bi{t). As soon as some a/(t) is found to be zero up to the 
required precision necessary to obtain the term of C/j, 
all ap{t) and bp{t), p > I become unecessary. 

The above procedure generalizes what is done in weak- 
coupling theory. There, Wick's theorem allows a resum- 
mation of one-particle reducible diagrams, which gives 
Dyson's equation. If the self-energy is LR, {i.e., has an 
underlying continued fraction structure) , the Green func- 
tion inherits this property due to the form of the weak- 
coupling free propagator. 

We were able to deduce from Eq. (||) the following 
continued fraction 



Qj{iU}n 
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2te(k)- _ 3/3*3 tanh (^^jc(k)/i 



Qdt^ 



iujn — 2*c(k)/(i— iujn + tc(k)/(i 



(3) 



, (2) 



which verifies the conditions CO, and has exactly the 
same Taylor expansion as Q up to order included. This 
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means that all the moments of Gj are the same as 
those of the exact solution except for terms of order i"*. 
Furthermore, any LR rational function sharing this prop- 
erty reduces to a continued fraction whose coefficients 
differ from those of Eq. only by terms smaller than 
the precision achieved here . 

Expansion to order for the half-filled HM suf- 
fices to exhibit both the Mott transition and the ef- 
fect of AF correlations on the spectral weight A(k, uj) = 
lim^_,o+ ~2 Im C/(k, a; -|- irj). There is no rigorous def- 
inition of the Mott transition in terms of one-particle 
properties, but one can use, as a heuristic criterion, the 
appearance of spectral weight at zero frequency. In the 
density of states N{uj) = J^^ A{k,u;)d'^k/{2iTY, as t in- 
creases from zero, the two symmetric Hubbard bands lo- 
cated at u and — u in the atomic limit widen, and even- 
tually mix for t beyond some critical value. The latter 
may be obtained by demanding that a pole of Q crosses 
the Fermi level for some k. For T ^ oo, the critical 
value of t is tc ^ u\/lTVTTT2¥ / {2dV3) ||. This 
gives Uc ~ 3.2i for (i = 1, to be compared with Uc — 3.5i 
found in the Hubbard-III |^ approximation. At finite T, 
we cannot calculate tc analytically, but Fig. sketches a 
numerical evaluation (for d = 1) in the (T, t) plane of the 
line where the gap vanishes. The value of grows upon 
lowering T, and there is no Mott transition at zero tem- 
perature, in agreement with the exact result of Ref . Q . 

The effects of AF correlations show up at low T, as 
illustrated in Fig. || by the plot of A{k,uj) for point B 
of Fig. |l| (k becomes k because we discuss the ID case 
for definiteness) . A{k,uj) has four delta peaks (a finite 
width T] is added for clarity) given by dispersion relations 
u>i{k), {i—l to 4 as in Fig. Q). The spectral weight is an 
even function of fc, and particle- hole symmetry ensures 
that A{k + TT, —Lo) = A{k,u}). While at small t and high 
T (point A), (fc) was minimum for fc = 0, when T is 
lowered down to point B, the minimum of W2(fc) moves 
continuously from k = towards fc = 7r/2 (Fig. ||), and 
peak 2 loses weight for values of k much smaller than 
tt/2. These changes reflect the AF short-range order that 
gradually builds up when T becomes smaller than the 
AF superexchange J — 2t^ /u of the equivalent t — J 
model. The approximate cell doubling in direct space 
translates into a nearly 7r-periodic dispersion for peak 2, 
although the 27r-periodicity of its weight and of uji{k) 
reminds us that the state remains paramagnetic. This is 
why we chose to define the AF crossover line of Fig. |l| 
as the points where fc = ceases to be the minimum of 
0^2 (fc)- In this regime, the width of band 2 is of order 
J = 21^ /u whatever the value of t, supporting the above 
interpretation. 

If we decrease T further from point B, we enter a 
regime that is beyond the domain of validity of our 
approach. Indeed, in contradiction with the results of 
Refs. the spectral function becomes similar to 
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FIG. 2. (above) Spectral function A{k, uj) for point B 
{t = 0.25w, T = 0.06« or U = 8t,T = 0.24t). rj (see text) was 
set to 0.02. (below) Dispersion relation of the peaks. 

that of free particles, following a —2tcosk dispersion, 
except for a gap at the Fermi energy for k ~ 7r/2. We ex- 
pect our expansion to be valid if the bi (t) 's in Eq. (^) are 
small compared to to, whose lowest-order value is u. This 
leads us to the conditions t/u < 0.5 and 3{t/u)^ < T/u, 
fulfilled by the points under the dashed line in Fig. 
However, these conditions may be too stringent because 
the t oo limit also happens to be correctly given by our 
solution Eq. (^). Furthermore, a free particle dispersion 
relation, with a gap opened at the Fermi level, is what is 
expected at large t and small T for an itinerant antiferro- 
magnet. Fig. g| (point C) illustrates this behaviour. The 
parameters have the same value as in the Monte-Carlo 
(MC) calculations of Ref. [| {U = At, l3 = 20/t). The 
general distribution of the spectral weight, and the dis- 
persion relation of the peaks Q are well accounted for 
by our solution. We believe that peak 1 contributes to 
the large uncertainty (due mainly to the Maximum En- 
tropy Method itself) on the maxima of A{k,uj) reported 
m Fig. 2 of Ref. (D for k near and n. For other values 
of fc, peak 1 could not be resolved in Ref. |^ because 
of its small weight and because of the magnitude of the 
time slice, unlikely to detect high-energy features. Thus, 
our results for point C appear correct. Our method defi- 
nitely fails in the shaded area of Fig. ^ where spin-charge 
separation occurs, but outside this region our solution is 
reliable under the dashed line, and uncontrolled (but not 
necessarily bad) above it. 



3 



-6 



CO (xO 

4 6 



7t/4- 



71/2 



371/4 









41 9-4-1- 












.A 








^ 


<y — } 




w — 


^ V nt 












-3 -2 -1 1 2 , ,3 

CO (X m) 






fll 




[2] 




f3l 




[4] 




[5] 




[6] 




[7] 


s 


[8] 


3^ 
X 


[9] 


2^ 




1 


[10] 




[11] 





[12] 


-1 


[13] 


-2 




-3 


[14] 




[15] 



FIG. 3. (above) Spectral function A{k, uj) for point C 
{t = 0.5u, T = 0.025m ot U = 4t, T = 0.05t, (3 = 20/t) with 
77 = 0.02. (below) Dispersion relation of the peaks. 



The cuprate chain material SrCu02 studied in Ref. ||] 
happens to fall in the shaded regime. Nevertheless, our 
results allow us to predict that features dispersing on a 
scale J, like peak 3 in Fig. should appear at 7r/2 < 
k < TT upon raising T. Hints of this finite T effect have 
already been seen in the "question-mark" features in Fig. 
1 of Ref. |]. 

In summary, we presented a general method for con- 
structing strong-coupling expansions and applied it to 
the half- filled Hubbard model. We showed how the Mott 
transition and AF correlations manifest themselves in the 
single-particle properties. Finally, we gained further in- 
sight into ongoing ARPES experiments on the propaga- 
tion of one hole in an AF correlated Mott insulator. Dop- 
ing and two-particle correlations are accessible within the 
same approach. 
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